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Abstract 

We propose a unified space-time picture of baryon stopping and quark-gluon 

plasma creation in ultrarelativistic heavy-ion collisions. It is assumed that 

the highly Lorentz contracted nuclei are decelerated by the coherent color 

field which is formed between them after they pass through each other. This 

process continues until the field is neutralized by the Schwinger mechanism. 

Conservation of energy and momentum allow us to calculate the energy losses 

of the nuclear slabs and the initial energy density of the quark-gluon plasma. 
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Ultrarelativistic heavy-ion collisions offer a unique possibility to study collective dynam- 
ical effects at the partonic level. Before a collision the partons are confined in coherent 
field configurations and the Lorentz contracted nuclei can propagate in the physical vacuum 
without distortion. After the nuclei collide, thousands of partons are quickly liberated. The 
nuclei are transformed into two partonic sheets which recede from each other, leaving behind 
the perturbative QCD vacuum and strong gluon fields [l] . Since color charges on the sheets 
are distributed stochastically, they generate chromoelectric fields which are nonuniform in 
the transverse plane. This field configuration may be envisaged as a collection of densely 
packed color flux tubes or strings stretched between the sheets. The energy accumulated 
in the coherent field is taken from the kinetic energy of the receding partonic sheets that 
causes their deceleration. At later times the coherent fields are neutralized via the Schwinger 
pair-production mechanism 0. 

In this letter we formulate a simple model to describe the collective deceleration of 
ultrarelativistic nuclei by the coherent color field. This model is similar in spirit to that 
proposed recently in ref. ||. Some geometrical and kinematical aspects of the model are 
also similar to those introduced in ref. |IJ]. All calculations below are performed in the 
center-of-rapidity frame where the colliding nuclei have initial rapidities of ±t/ . 

We consider only beam energies so great that the nuclei can be thought of as very thin, 
Lorentz contracted, sheets. Each sheet is divided into many small elements, "slabs" , of unit 
transverse area labeled by an index a where a = p for the projectile nucleus and a = t for 
the target nucleus. Each slab is characterized by a baryon number N a which is assumed to 
be strictly conserved. This number can be expressed through the slab thickness, l a , in the 
rest frame of the respective nucleus, 



where p a (b, z) is the density distribution in nucleus a, b is the position of the slab in the 
impact parameter plane, and po ~ 0.15 fm -3 is the equilibrium nuclear density. For a 
uniform density distribution 1(b) = 2\/R 2 — b 2 , where R = r A 1//3 is the nuclear radius and 
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(1) 



r = (47rp /3) -1 ^ 3 ~ 1.17 fm. We assume that before and after nuclear overlap, which takes 
place at t — 0, each slab propagates as a rigid body along the beam axes z. The energy 
and momentum of slab a is parameterized in terms of its proper energy per baryon e a and 
longitudinal rapidity y a (for details see ref. M), 

E a = N a e a cosh y a , P a = N a e a sinh y a . (2) 

In general, e a may differ from the energy per baryon in normal nuclei, Eq ~ m^r, the nucleon 
mass, due to hard parton production at t — 0. 

As a result of soft parton-parton interactions the slabs acquire color charges which gen- 
erate the chromoelectric field or several strings in the space between them, see Fig. 1(a). 
For simplicity we disregard here a short time delay which is needed for the formation of the 
coherent field. At t > the trajectories of the projectile and target slabs, z p (t) and z t (t), 
are affected by the energy and momentum losses for stretching the strings. If the string 
tension (energy per unit length) for individual strings is a, and the number of strings per 
unit transverse area is n, then the potential energy stored in strings is V(z) = na\z p — z t \ 
(for definiteness we assume that z p > z t ). Accordingly the force or, in the present context, 
the longitudinal pressure exerted by strings on slab a is —dV/dz a = =Fna, and therefore 
Newton's equation of motion has the form 

dP n 



dt 



Tna . (3) 



From here on the upper and lower signs correspond to a = p and a = t, respectively. The 
net force acting on both slabs is of course zero and the total energy of the slabs+field system 
is conserved. Realizing that the energy lost by a slab while traversing distance dz a is equal 
to ^nadz a , one can also write the energy conservation equation as 

dE a 

These equations are identical to those describing the motion of massive capacitor plates 
due to the action of a uniform (chromo)electric field of a capacitor. This is not surprising 



because the action of many collinear flux tubes, or strings, is equivalent to the action of a 
uniform field E with the same energy density, 6f — E 2 /2 = no . It is important to note that 
Sf is a Lorentz invariant quantity. 

Using eqs. (||) and (f|) and the definition of slab velocity, 

dz P 

one finds that P a dP a = E a dE a . This means that the internal energies of the slabs do not 
change in the course of deceleration, that is E% — P% = N%£% = constant. The solution to 
equations @ and (f|) for the initial condition y a (0) = ±yo, z a (0) = are easily found to be: 

sinh y a = ± ^sinh y - > (6) 

cosh y a = cosh y =p , (7) 

where A a = (e a p /ef) l a is the characteristic deceleration length. 

Eliminating y a from these equations one obtains the slab trajectories 

^coshyo T y) ~ ( sinn ?/o - y) = 1 • ( 8 ) 

These are parts of hyperbolae shown in Fig. 2. Initially the slab trajectories are very close to 
the light cone but later on they increasingly deviate from it. If nothing else were to happen, 
the slabs would reach their respective turning points at t = X a sinh?/ , z = ±X a (coshyo — !)> 
and then reverse direction. As an example, Fig. 2(b) shoes the projectile trajectory with a 
turning point. This is the "yo-yo" type of motion well known in string models. However, such 
dynamics is very unlikely in nuclear collisions because of irreversible processes associated 
with string decay. 

When the strings become long enough quark- ant iquark and gluon pairs can be produced 
from vacuum via the Schwinger mechanism [Q]. Their color charges will screen the chromo- 
electric fields and eventually neutralize them. For various idealized situations these processes 
were studied intensively during the last decade 0]. Here we want to avoid numerical sim- 
ulations of the complicated plasma- field dynamics, which is not very well understood yet. 



Instead we adopt a simplified picture assuming that the strings decay suddenly at a certain 
proper time r = r ~ 1 fm, where r = \Jt 2 — z 2 . As a result, the partonic plasma is created 
at the hyperbola t = Tq (see Fig. 2). As we will see below, even this simplified picture 
allows us to make nontrivial predictions concerning the energy density and velocity field of 
the plasma at r = tq. 

As shown in Fig. 1(b), the structure of the system changes dramatically at t > Tq. Now 
the chromoelectric fields are absent in the middle but still remain in the regions adjacent to 
the projectile and target slabs. The gap between these two regions is filled by the partonic 
plasma. As time progresses the boundaries between the plasma and coherent fields move 
closer to the slabs, and the gap expands. Since in a time interval dt the boundaries traverse a 
distance dz along the hyperbola y/t 2 — z 2 = r , their velocities are dz/dt = tj ' \jt 2 — Tq > 1. 
Although this velocity is superluminal, it has no physical significance since the boundary 
is not a material object. The plasma itself moves in accordance with the laws of special 
relativity. To see this, let us write the energy-momentum tensor of the plasma in the 
standard form 

= (e + p)uV - pg^ , (9) 

where e and p are the proper energy density and pressure, is the collective 4-velocity 
of the plasma, and g M ^=diag(l,-l,-l,-l). At its creation the plasma is assumed to have 
negligible transverse velocity; collective transverse expansion develops at later times. Thus 
the 4-velocity can be chosen as = 7(1, 0, 0, v) where v is the longitudinal velocity and 
7 = — v 2 . Now we can write equations expressing the conservation of energy and 

momentum across the boundary between the plasma and the chromoelectric field. In time 
dt, when the boundary moves a distance dz, the energy per unit area subtracted from the 
field is €fdz. This change must be equal to the energy of a newly produced plasma slice dz: 
T 00 dz = €fdz. Using eq. @ this gives 

(e + p) 1 2 -p = e f . (10) 
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Since the field exerts a force na = e/ on the plasma charges, the change of momentum in 
the plasma slice dz is T 03 dz = £fdt. This leads to 

(e + p)-f 2 v = e f -t = e f Z - . (11) 

In the last equality we have used the condition that the boundary moves along the hyperbola 



\/t 2 — z 2 = r so that tdt = zdz. The two equations ( |TDD and (|TT|) allow us to find two 
quantities, e and v , characterizing the plasma at r = tq. Especially simple results are 
obtained in the case of free streaming, p = (dust equation of state), which seems most 
appropriate for the early stages of the plasma evolution. By dividing eq. (|TTD by eq. fliCf ) 

we get 

v(Vt 2 - z 2 = t ) = - = tan 77 , (12) 

t 

e(V^^=r ) = % = ^- , (13) 
7 cosh 77 

where the space-time (pseudo)rapidity 77 = | In ffzf) has been introduced. The first formula 
gives exactly the velocity field postulated in Bjorken's scaling hydrodynamics 0, but here 
it follows from the conservation laws at the plasma boundary. The second equation shows 
that it is the global energy- density T 00 = e7 2 and not the rest-frame energy density e which 
should be constant at the hypersurface r = tq where the plasma is created. As follows from 
eq. fljlf ) the actual proper energy density at r = r decreases rather rapidly when going 
away from the central point 77 = 0. These predictions should be used as the initial conditions 
for further kinetic or hydrodynamical simulations of the plasma evolution. 

The trajectories of the projectile and target slabs are not affected by the string decay until 
they intersect the hyperbola t = tq. It is clear from the previous discussion that the region 
occupied by the plasma expands faster than the receding slabs. Thus the plasma eventually 
eats up the strings at time t* when the slab trajectories intersect with the hyperbola r = r , 
as shown in Fig. 2. At this time the chromoelectric fields are fully neutralized and no new 
plasma is produced any more (see Fig. 1(c)). Substituting z* = z a (t*) in eq. ([?[) we obtain 
the final gamma-factors and rapidities of the slabs, 



7a = coshy* == 7o 



1 r ° I 

1 - T- \ V 0< 



TO 



(14) 



A a ^ u \| AXl 2X a 

where 7 = coshyo and t> = tanhy . We see that 7* is entirely determined by a single 
combination of parameters, r /X a = (efTo)/(e a p l a ). Since e a ~ m^y is not expected to vary 
much in the course of the reaction and l a is given by the geometry, the combination e/To is 
the only essential parameter of the model which determines the rapidity loss by the nuclei. 
It is quite natural that this is the same quantity which gives the initial transverse energy of 
the plasma at central pseudorapidity 77 = 0. As follows from eq. (|i~3l), the energy density 
of the plasma at rj = and r = r is equal to ej. The total energy of the plasma in a slice 
dz = Todrj around rj = is obtained by integrating e/(b) over transverse area of the reaction 
zone. This gives for central collisions 

^ = 0) = rf( e/ )r Q , (15) 

where (ej) is the energy density of the chromoelectric field averaged over the transverse area. 

As our estimates show, it is most likely that A a r . In this case the partonic slabs 
are significantly moderated at the time t* a . Indeed, in this limit t* a ~ z* ~ A a 7o, which is 
close to the turning point. From eq. (0) the condition of complete stopping (7* = 1) is 



r o/A a = ^2(70 — 1). In collisions of asymmetric slabs the shorter slab may even reverse its 
motion in the direction of the longer slab, as illustrated in Fig. 2(b). Therefore, the proposed 
scenario provides an efficient mechanism for baryon stopping in ultrarelativistic heavy-ion 
collisions^. The baryon rapidity distribution is obtained by summing up contributions of all 
pairs of slabs and integrating over impact parameters. Our prediction is that in the first 
approximation this will be a superposition of two peaks centered at rapidities y* and y\ 



given by eq. (14). In a more refined approach one should also study what happens to the 



baryonic slabs after they are hit by the plasma. It is clear that the plasma wind will cause 



1 One should keep in mind, however, that even in the case of complete stopping the projectile 
and target slabs will be separated in space by a distance of about (X p + A^o- 
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drift and diffusion of the baryon charge in rapidity space. As a result the projectile and 
target peaks will be shifted towards initial rapidities and smeared out. 

The process of plasma creation described above is similar to the breakdown of a dielectric 
in a strong electric field. The electromagnetic plasma produced in this case shows up as a 
spark or lightning. Since here we are dealing with the QCD plasma, this phenomenon can 
be referred to as "QCD lightning". 

A direct calculation of the energy density e/ accumulated in the chromoelectric field at 
the early stage of a heavy-ion collision is problematic at present. Therefore, for our estimates 
we use a simple parameterization which is motivated by several model calculations. 



Here the second factor is responsible for the incident energy dependence, ^/s is the cm. 
collision energy per nucleon pair. The exponent a ~ 0.3 follows from the low-x behavior 
of the nuclear structure function or parton density J7|. The third factor accounts for the 
geometry of the overlap zone. The exponent f3 relates the number of strings produced to the 
number of binary parton-parton collisions, which is proportional to N p N t . For convenience 
we have normalized N p and N t by the mean areal baryon density in the proton, N « 0.5 
fm -2 . In the case of uncorrelated strings, as in elementary pp or pp collisions (5 ~ 1, while in 
the case of strong string overlap or percolation (3 ~ 0.5 ||. For the estimates below we use 
the value j3 = 0.5, also used in ref. ||, which seems appropriate for heavy- ion collisions at 
RHIC and LHC energies. It is assumed that this parameterization applies above a certain 
cm. energy squared, s > Sq. The parameter eo can be determined at s = Sq and then used 
for higher energies. With this choice we get 



where a = t for a = p and vice versa. This formula shows that in an asymmetric slab- 
slab collision (l a < /„) the smaller slab will decelerate faster than the bigger one. This is 
a natural result since according to eq. (|]) their energy losses are equal. For equal slabs 




(16) 




(17) 



S 



(Ip = h = I) the deceleration constant is independent of their lengths. This means that in 
a central collision the whole partonic sheet will decelerate uniformly. On the other hand, 
according to eq. (|16"1), the energy density of the chromoelectric field decreases when going 
from the symmetry axis to the periphery, e/(b) oc /(b). 

We shall now discuss phenomenological implications of our model in light of recent 
RHIC data for central Au+Au collisions at y/s = 130 GeV (yo = 4.94). As found by 
the PHENIX collaboration |§, the transverse energy in the central pseudorapidity window 
is cIEt / dr/l^o m 600 GeV. For the free streaming plasma (p = 0) e(r)r = e(r )r , and 
dEx I cZ?7 1 T7=o is independent of r. Then from eq. (|i~5"D with irR 2 = 148 fm 2 we immediately 
obtain (e/)r ~ 3.9 GeV/fm 2 . If the pressure effects are fully included, p = e/3, the trans- 
verse energy drops as r _1//3 , and in order to obtain the observed dE^jdr] one should increase 
the initial value by a factor (tj/tq) 1 / 3 ~ 2, where r and tj are the initial and final time of 
hydrodynamic expansion. In our estimates we take for (e/)ro the value 6 GeV/fm 2 which 
is in between these two extremes. This value is in qualitative agreement with results of 
microscopic simulations within the Quark- Gluon- String Model [fL0| . 

Now we can estimate the final rapidities of Au nuclei by considering the collision of 
representative slabs with lengths l p = l t = (I) = 4R/3 averaged over the transverse plane. 
Below we omit index a for simplicity. The proper energy per baryon in partonic slabs can be 
approximated by rriT = \fm% + (pt) 2 , where (px) ~ 0.5 GeV/c as measured by the STAR 
Collaboration [|ll|] . To ensure energy conservation at t = we have reduced 70 by the factor 
in n j uij'. With these inputs we find from eq. (|17]) that To/ A = 4.1 or the characteristic 
deceleration length A = 0.24r 0. From eq. ([14]) we get the final projectile and target 
rapidities, y* ~ ±1.9. This corresponds to the final cm. energy of about 3.5 GeV per 
baryon compared to the initial energy of 65 GeV per nucleon. The difference is transferred 



2 It is interesting to note that with to = 1 fm this would correspond to the total deceleration time 
and distance t* ~ z* ~ 14 fm. 
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into the quark-gluon plasma. This is a tremendous energy loss which can explain the rather 
high degree of baryon stopping observed at RHIC. Unfortunately, the net-baryon rapidity 
distributions for this reaction are not available yet. Preliminary BRAHMS data show that 
the net-proton rapidity distribution has a dip at y — and two peaks at y — ±(2 4- 3) flT2|| . 

Finally we use eq. fllTf ) to make predictions for the full RHIC energy ^/s = 200 GeV. 
Extrapolating from ^/s = 130 GeV to ^/s = 200 GeV, we get (e/)r = 6.8 GeV/fm 2 and 
y* = ±2.1. If instead we would use the values e(ro) = 66.5 GeV/fm 3 and tq = 0.3 fm 



given in ref. pM , we would predict that the baryonic slabs reach the turning point and 
reverse the direction of motion. Strict correlation between the transverse energy produced 
at central rapidities and the degree of baryon stopping is the most important prediction of 
our approach. The predicted strong deceleration of the baryon as well as the electric charge 
of ultrarelativistic nuclei makes it very promising to search for photon bremsstrahlung at 



collider energies (see ref. |14]] and references therein). 

In conclusion, the collective dynamics of partonic slabs has been studied within a simple 
model incorporating strong chromoelectric fields generated early in the reaction. The model 
has two essential physical parameters: the energy density accumulated in the chromoelectric 
field, ef, and the proper time required for its neutralization, r . Applying conservation of 
energy and momentum, we have established a direct correspondence between the rapidity 
shifts of the projectile and target nuclei and the transverse energy of the produced quark- 
gluon plasma. Using data from recent RHIC experiments we have estimated the energy 
density accumulated in the chromoelectric field and the resulting rapidity shift in the net 
baryon distribution. 
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edges the kind hospitality of the Nuclear Theory Group, University of Minnesota. This 
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FIGURES 
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FIG. 1. Spatial view of a slab-slab collision at different times: (a) < t < To, (b) to < t < t* 
and (c) t >t*. The black boxes represent the receding projectile and target slabs. The horizontal 
straight lines indicate strings. The dotted areas show the regions occupied by the quark-gluon 
plasma. 
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FIG. 2. Schematic space-time representation of a symmetric (a) and asymmetric (b) slab-slab 
collision in the t — z plane. The projectile and target slab trajectories are shown by thick solid 
lines which start at the origin and terminate at points P and T, respectively. The quark-gluon 
plasma is produced at the portion of the hyperbola r = tq between these points (thick solid line). 
Horizontal solid lines represent strings. For symmetric collisions t* = t^=t*,z* = —z% = z*. 
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